



double v_sprod( double *a, double *b)
/* *********************************************************   */
/* * Skalarprodukt: S=A*B         alfa = acos(S/(|A|*|B|)) *   */
/* *********************************************************   */
{
  return(a[0]*b[0]+a[1]*b[1]+a[2]*b[2]);
}
